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The ordering of the Ising model on a pyrochlore lattice interacting via the long-range 
RKKY interaction, which models a metallic pyrochlore magnet such as Pr 2 Ir 2 07, is stud- 
ied by Monte Carlo simulations. Depending on the parameter kp representing the Fermi 
wavevector, the model exhibits rich ordering behaviors. 
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Recently, the ordering of pyrochlore magnets consisting of corner-sharing network of tetra- 
hedra has attracted much attention due to an intrinsic effect of geometrical frustration. Inten- 
sive studies were performed in particular on "spin-ice" magnets, e.g., Ho2Ti207 and Dy2Ti207, 
which possesses a strong Ising- like < 111 > magnetic anisotropy and exhibits an intriguing 
glassy behavior analogous to the one associated with the proton ordering of water ice. 1-3 ) 
In spin-ice magnets, the local spin configuration at each tetrahedron is known to be the so- 
called "2-in 2-out" configuration, which is stabilized by an effective ferromagnetic interaction. 
Indeed, this "2-in 2-out" condition is reminiscent of the "ice-rule" condition familiar in the 
proton ordering of water ice. While an earlier analysis on spin ice was made based on the 
nearest-neighbor ferromagnetic exchange model, 1 ) the dominant interaction in real spin-ice 
magnets is the long-range dipolar interaction, which falls off as 1/r 3 with distance r. 3 ~ 6 ) 
Thanks to extensive experimental and theoretical studies, the experimental results on spin ice 
magnets now seem to be relatively well understood on the basis of the dipolar spin-ice model, 
i.e., the pyrochlore Ising model with the long-range dipolar interaction. 

Typical spin-ice magnets Ho2Ti207 and Dy2Ti2C>7 are insulators. Ordering properties of 
the corresponding metallic pyrochlore magnets are then of special interest. 3 ) Recently, Nakat- 
suji and collaborators reported that a metallic Ising-like pyrochlore magnet P^I^Oy exhibits 
an interesting behavior at low temperatures. 7 ) The magnetism of this compound is borne by 
the localized moment of Pr 3+ which possesses a strong < 111 > Ising-like magnetic anisotropy. 
Reflecting the metallic nature of the compound, the dominant magnetic interaction is expected 
to be the long-range RKKY interaction mediated by the Ir 4+ conduction electrons, rather than 
the dipolar interaction. Although the Curie- Weiss constant of P^L^Oy determined from the 
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high-temperature susceptibility measurement is antiferromagnetic #cw — — 20K, and the an- 
tiferromagnetic interaction usually prefers "all-in all-out" structure accompanying no frustra- 
tion, the short-range magnetic order of this compound realized at lower temperatures appears 
to be "2-in 2-out" similarly to that of spin ice, characteristic of the effective ferromagnetic 
interaction. 

While in spin-ice magnets the observed ferromagnetic Curie- Weiss temperature is entirely 
consistent with the "2-in 2-out" local spin structure, 1 ) it remain a bit puzzling in P^h^Oy 
how the "2-in 2-out" structure, which is usually an attribute of the ferromagnetic effective 
interaction, is realized under the antiferromagnetic ©cw- On decreasing the temperature, 
Pr2lr207 exhibits a spin- liquid- like behavior accompanied by the Kondo-like effect observed 
in its resistivity, until it eventually exhibits a weak spin freezing at a low temperature T ~ 
0.12K. 7 ) 

Under such circumstances, to better understand the ordering of metallic pyrochlore mag- 
nets, it would be desirable to study the nature of the magnetic ordering of the pyrochlore Ising 
model interacting via the long-range RKKY interaction, and compare it with the property of 
the dipolar spin-ice model. 

We consider the Hamiltonian 



where Si denotes a three-component unit vector along the < 111 > direction located at the site 
i of the pyrochlore lattice, pointing either parallel or antiparallel with the vector connecting 
the site i and the center of a tetrahedron. The coupling between the spins at the sites i and 
j is the long-range RKKY interaction given by 



where is the distance between the sites i and j, kp the Fermi wavevector, and the sum is 
taken over all spin pairs on the pyrochlore lattice. 

The choice of the standard RKKY interaction (2) is justified when conduction electrons 
are sufficiently extended and are coupled with localized moments via the on-site exchange 
interaction. It is not clear at the present stage to what extent this assumption is satisfied in 
real metallic pyrochlore magnets such as P^fr^Oy. If, e.g., the Ir 4+ conduction electron has 
negligible amplitude at the Pr 3+ site, the d-f superexchange might play a role modifying the 
standard form of the RKKY interaction (2) , although even in such a situation the qualitative 
feature of Eq.(2), i.e., the long-range and the oscillating nature of the effective interaction, is 
still expected to hold. Furthermore, in understanding certain properties of real systems, the 
itinerant character of electrons, not captured by the spin model (1), might become important. 
In the following, with these possible limitations of the model in mind, we study the order- 
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ing properties of the RKKY Ising model (l)-(2), as a first step to understand the ordering 
properties of metallic pyrochlore magnets. 

The model contains one parameter, the Fermi wavevector In case of P^I^Oy, ex- 
perimentally determined carrier density n ~ 4.13 x 10 21 cm~ 3 and the lattice constant 
a ~ 10.39 x 10~ 8 cm 7 ) yield the dimensionless Fermi wavevector kp — 27r/1.218 in units 
of a -1 . Our length unit is taken to be the spacing of the cubic unit cell a. The lattice con- 
tains N = 16L 3 spins (the cubic unit cell contains 16 sites). Periodic boundary conditions are 
applied. 

In order to study the ordering properties of the model, we perform Monte Carlo simulations 
on finite lattices, with L in the range of 2 < L < 6. In order to take account of the long- 
range nature of the RKKY interaction properly, we employ the Ewald sum technique, i.e., 
the RKKY interaction is summed over to infinity under the condition that the images of the 
original lattice are repeated periodically in all directions. 8 ) Typically, each run consists of 
5 x 10 4 Monte Carlo sweeps, several independent runs being made to guarantee the statistical 
accuracy of the data. Many of the data given below are taken at the "experimental" &f- value, 
kp = 27r/1.218, while other /cp-values are also investigated. 

First, we concentrate on the "experimental" case of Hf = 27r/1.218. In Fig. 1(a), we show 
the temperature dependence of the specific heat per spin. The data exhibit a sharp peak at 
around T = T c ~ 16 (in units of Jo). The peak grows rapidly with the size L, suggesting 
that the transition is of first-order. Indeed, as shown in the inset, the energy histogram, i.e., 
the distribution function of the energy per spin e, exhibits a double-peak structure around 
T c , which grows rapidly as L increases. In the inset of Fig. 1(b), the inverse susceptibility 
X 1 = (Xxx + Xyy + Xzz)" 1 i s shown as a function of temperature. The ordered state has no 
net magnetization. The estimated Curie- Weiss constant is negative, @cw — — 32, apparently 
suggesting the occurrence of an antiferromagnetic effective interaction. 

Naively, one may think that the antiferromagnetic Ocw is associated with the "all-in 
all-out" local spin structure. In fact, this is not the case here. The inset of Fig. 1(b) exhibits 
the temperature dependence of the ratios of the "2-in 2-out" structure, the "3(l)-in 1(3)- 
out" structure and the "all-in all-out" structure, which represent the probability that each 
type of local spin structure is realized at each tetrahedron. As can be seen from the figure, 
the local spin structure below T c is "2-in 2-out" , not "all-in all-out" . The appearance of the 
ferromagnetic "2-in 2-out" spin structure under the antiferromagnetic 6cw is consistent with 
the experimental observation for P^I^Oy. 7 ) In fact, as we shall see below, such an apparently 
contrasting behavior can be understood as a consequence of the long-range and oscillating 
nature of the RKKY interaction. 

To further clarify the nature of the ordered state, we calculate the spin structure factor 
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Fig. 1. (Color online) The temperature dependence of the specific heat per spin (a), and of the 
temperature dependence of the ratio r of the "2-in 2-out", "3(l)-in l(3)-out" and "all-in all-out" 
local spin structures at each tetrahedron (b), for the parameter kp — 27r/1.218. In the inset of (a), 
the energy histogram is shown at a temperature close to the transition temperature, at T = 15.750 
(L = 4) and at T = 15.875 (L = 6). The inset of (b) exhibits the inverse susceptibility per spin 
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F(q), i.e., the thermal average of the squared Fourier amplitude defined by 

i? ^ = ^<IE^ e ^l 2 >' (3) 

j 

where the wavevector is given by q = —{h, k, I) with —2 < h, k, I < 2. In Fig. 2(a), we show the 
structure factor in the (h, h, I) plane calculated at a low temperature T = 2, which has been 
averaged over all equivalent directions of (±/i, ±/i, ±Z), (±/, ±/i, ±/i) and (±/i, ±Z, ±h). As can 
be seen from the figure, the spin structure factor exhibits a sharp Bragg peak at (0,0,1), 
or at the equivalent ones, indicating that the magnetic long-range order sets in below T c , 
the associated spin order being characterized by the wavevector (0,0, ^l). We have confirmed 
that the Bragg peaks appear at these spots irrespective of L, and sharpen with L. It is not a 
finite-size effect. 
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Fig. 2. (Color online) The spin structure factor F(q) in the (h, I) plane at a temperature T — 2 for 
kp = 27r/1.218 (a), and for kp = 27r/1.13 (b). The system size is L = 6. 



In fact, such an ordered state is the same as the one observed in the dipolar spin-ice 
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model in its equilibrium. 9 ) In the case of the dipolar spin-ice model, the (0,0,^) ordered 
state is observed only in a special type of simulation devised to promote equilibration (loop 
algorithm), not in the standard single spin-flip simulation: 6 ) In the latter case, the system 
gets out of equilibrium leading to the glassy, spin-ice behavior. In our present RKKY model, 
by contrast, the (0,0, ^) long-range order is achieved more readily, even without using the 
loop algorithm. 

The (0,0,—) ordered state is six-fold degenerate reflecting the cubic symmetry of the 
lattice. The order parameter possesses a structure equivalent to the one of the ferromagnetic 
cubic model. Since the latter model is known to exhibit a first-order transition in three di- 
mensions, 10 ) the observed first-order nature of the transition is naturally understood as a 
consequence of the occurrence of the (0,0, ordered state below T c . 

In order to examine the stability of the ordering against the possible change of the param- 
eter kp, we also simulate the model with other values of kp. For the wavevectors kp = 27r/1.15 
and 27r/1.4 close to the case studied above, we have observed essentially the same ordering 
behavior, i.e., a first-order transition into the (0,0, — ) ordered state consisting of the "2-in 
2-out" local spin structure, which is accompanied by the negative 0cw- 

If one further increases or decreases the A; jr- value, somewhat different behaviors arise. 
With increasing the Ap-value, |©cw| decreases its magnitude, eventually changing its sign 
from negative to positive at around kp ~ 27r/1.14. At kp = 27r/1.13, Ocw becomes positive 
(ferromagnetic), while the model exhibits a first-order transition into the long-range ordered 
state with keeping the "2-in 2-out" local spin structure. The ordered state is characterized 
by the wavevector (-,-,-) which is distinct from the (0,0,—) ordered state observed at 
kp = 27r/1.218. The computed spin structure factor at kp = 27r/1.13 is shown in Fig. 2(b). 

On the other hand, with decreasing the fcF-value from kp = 27r/1.218, the sign of 0cw 
is kept to be negative and the local spin structure remains to be "2-in 2-out", while at 
kp = 2n/1.7 the model exhibits a first-order transition into the more complex long-range 
ordered state characterized by the wavevectors with all three spin components different from 
each other, h^k^l. 

In Fig.3, we summarize the ordering behavior of the model for a wider range of kp. The 
sign of the Curie- Weiss constant 6cw determined from the computed \ l an d the type of 
the local spin structure, i.e., either ferromagnetic "2-in 2-out" or antiferromagnetic "all-in 
all-out", are shown. We emphasize that, at any kp-value studied, we have never observed a 
spin-liquid-like behavior in which the spin long-range order is suppressed down to T = or 
temperatures much lower than 0cw- Further details will be reported elsewhere. 

The origin of the sign of the effective couplings can be understood as follows. First, the 
type of the local spin order, i.e., either "2-in 2-out" or "all-in all-out", turns out to be primarily 
dictated by the sign of the nearest-neighbor part of the RKKY interaction. To see this, we 
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include in Fig. 3 the sign of the nearest-neighbor part of the RKKY interaction as a function 
of kp, which turns out to reproduce the sign of the effective coupling associated with the local 
spin structure. 

By contrast, distant-neighbor interactions also contribute to its Curie- Weiss constant. This 
can be seen most easily from the lowest-order high-temperature expansion, which yields 

jGsub(i) j&ub(i) 

where the summation j € sub(i) is taken over all sites j belonging to the same fee sublattice 
as the site i, while the summation j € sub(z) taken over all sites j belonging to other three fee 
sublattices. The factor 1/9 in the second term arises due to the mutual canting of the < 111 > 
easy axes at different sublattices. The sign of eq.(4) is also shown in Fig. 3. As can be seen 
from the figure, it reproduces the sign of ©cw determined from our simulation reasonably 
well. Thus, in case of the long-range interaction like the RKKY interaction, the local spin 
structure at each tetrahedron is primarily dictated by its nearest-neighbor part, while the 
distant-neighbor interactions also contribute to the Curie- Weiss constant. 
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Fig. 3. (Color online) The A^-dependence of the sign of the Curie- Weiss constant Ocw and the 
sign of the effective coupling associated with the type of the local spin structure, i.e., either 
ferromagnetic "2-in 2-out" or antiferromagnetic "all-in all-out" . The information of the wavevector 
q characterizing the ordered state, whether it is q = or q ^ 0, is also given. For comparison, the 
sign of the nearest-neighbor part of the RKKY interaction J nn and the sign of Gcw calculated 
from the high-temperature expansion (H.T.ex.), eq.(4), are given. 

In the case of dipolar spin ice, the in-field properties turned out to be quite rich and 
interesting. 11 ^ In view of this, we also investigate the ordering behavior of our RKKY model 
under magnetic fields applied along the [111] direction for the "experimental" case of kF = 
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Fig. 4. (Color online) (a) The magnetization per spin plotted versus the [111] field intensity at a 
low temperature T = 0.5 for kp = 27r/1.218, together with the ratio r of the "3(l)-in l(3)-out" 
local spin structure. The lattice size is L — 4 and 6. The inset is a phase diagram in the field- 
temperature plane. In the shaded region, the detailed phase structure cannot be identified due 
to strong finite-size effects, (b) Spin pattern on a (111) kagome layer at the intermediate plateau 
state with m = 3/8. The black (white) circle represents the "in" ("out") spin state. The red or 
thick (blue or thin) triangle represents a face of the "3-in 1-out" ("2-in 2-out") tetrahedron. 

27r/1.218. The calculated magnetization per spin m is shown in Fig. 4(a). In addition to the 
low-field plateau at m = corresponding to the (0,0,^) ordered state and the high-field 
one at m = 1/2, there appears for intermediate fields another plateau at m = 3/8. Note that 
the pyrochlore lattice can be viewed as an alternate stack of a sparse triangular layer and a 
kagome layer. Under sufficiently strong [111] magnetic field, all tetrahedra take "3(l)-in 1(3)- 
out" configurations where all spins on a triangular layer are aligned parallel with the [111] 
field, with the saturation magnetization m = m s = 1/2. In the intermediate plateau state, a 
quarter of tetrahedra takes "3(l)-in l(3)-out" configuration, while other 3/4 take "2-in 2-out" 
configuration: See Fig. 4(a). These "3-in 1-out" tetrahedra are arranged periodically on the 
pyrochlore lattice, i.e., these tetrahedra form a 2 x 2 triangular superlattice on a (111) kagome 
layer as depicted in Fig. 4(b). In the inset of Fig. 4(a), we sketch a phase diagram of the model 
for kp = 27r/1.218 in the field-temperature plane. In the shaded region of the phase diagram, 
due to the possible the incommensurability effect, complicated behaviors accompanied by 
various narrow plateaus are observed, which, however, are strongly size-dependent in the 
range of small sizes studied here. 

Finally, on the basis of our present finding for the RKKY pyrochlore Ising model, we wish 
to discuss recent experiments on metallic pyrochlores. Our present result is consistent with 
the experimental result on P^I^Cv in that the local spin structure is "2-in 2-out" whereas 
the Curie- Weiss temperature is antiferromagnetic. 7 ^ By contrast, a first-order transition into 
the (0,0, ^) ordered state as observed in our simulation is not observed experimentally. If 
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one deduces the transition temperature by fitting the calculated 0cw to the experimental 
6cw, one gets an estimate of T c ~ 10K. Experimentally, below 10K, Pr2ir207 gets into the 
Kondo regime. This coincidence strongly suggests that the hybridization effect between the 
itinerant Ir conduction electrons and the localized Pr moments, not just the effect of mag- 
netic frustration, is essential in realizing the spin-liquid-like behavior observed in P^h^Oy at 
low temperatures. Meanwhile, small difference in samples may possibly stabilize the type of 
long-range ordered state revealed here. Indeed, a series of pyrochlore magnets, e.g., pyrochlore 
iridates i?2l r 20"7 including both metallic and insulating, exhibit a variety of ordering behav- 
iors, 12 ' 13 ^ and our present result on the RKKY pyrochlore Ising model might serve as a useful 
reference in interpreting the experimental data. 

The authors are thankful to Dr. S. Nakatsuji, Dr. Y. Machida, Dr. S. Onoda, Dr. M. 
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